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A method  for  obtaining  asymptotic  solutions  of  the  unstable 
resonator  integral  equation  which  is  valid  for  all  values  of 
the  magnification  was  developed.  Approximations  were  made 
on  the  Greens  functions  rather  than  the  eigenmodes,  leading 
to  results  which  are  easily  generalized  to  different  mirror 
geometries.  “Diffraction  dominated  eigenmodes*  for  reson- 
ators where  each  ray  escapes  after  a few  transits  were 
differentiated  from  ‘waveguide  dominated  eigenmodes*  which 
are  obtained  for  cavities  with  a large  number  of  transits 
per  ray.  The  solutions  obtained  were  seen  to  agree  in  the 
appropriate  limits  with  other  asymptotic  solutions,  numeri- 
cal results,  and  geometric  optics  predictions.  To  include 
the  effects  of  gain,  the  unstable  resonator  equation  was 
derived  from  Maxwell's  equations  in  a polarizable  medium. 

The  resulting  equations  have  the  same  structure  as  the  empty 
resonator  equation,  and  similar  approximations  can  be  used. 

Some  features  of  the  effects  of  saturation  on  the  eigenmodes 
of  an  unstable  resonator  were  considered^ 
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STATEMENT  OF  WORK 


contractor  shall  furnish  scientific  effort  during  the  period  and 
at  the  level  indica'T>d>4n^paragraph  la  of 'Section _H,  together  with  all  related 
services,  facilities,  suppl ies>«i^material s,  needed  to  conduct-the  following 
research. 


a.  Develop  and  apply  analytic  techniques  for  determining 
modes  and  eigenvalues  of  cylindriical  unstable  optical  resonators. 


b.  For  systems  which  display  azimuthal  symmetry,  develop  and  apply 
two  different  techniques:  (1)  an  asymptotic  differential  equation  approach 
and  (2)  an  approach  that  is  valid  for  situations  in  which  the  resonator 
magnification  is  large. 

c.  Develop  a technique  that  is  designed  to  treat  cylindrical 
resonators  in  which  the  system's  azimuthal  symmetry  has  been  destroyed  by 
optical  perturbations. 

d.  Develop  techniques  that  are  specifically  designed  for  resonators 
with  large  Fresnel  numbers. 

e.  Develop  a technique  to  incorporate  gain  in  unstable  resonator 


calculations. 


ANNUAL  REPORT 


Recently  there  has  been  considerable  interest  in  developing  analytic 
approaches  to  unstable  optical  resonators,  especially  devices  with  large 
Fresnel  numbers.  Two  reasons  for  this  interest  are: 

(1)  The  fundamental  physics  of  these  resonators  is  not  yet  fully 
* — understood. 

N 

' — - (2)  There  is  a need  for  rapid  and  accurate  techniques  for  calcula- 

ting the  eigenvalues  and  eigenmodes  of  large  Fresnel  number  un- 
stable resonators.  Such  devices  are  of  greatest  practical 
importance,  and  are  most  tractable  to  asymptotic  techniques.  We 
note  that  computer  generated  solutions  of  the  Fox-Li  resonator 
equation  become  more  costly  and  less  accurate  as  the  Fresnel 
number  increases. 

In  the  last  year  we  have  developed  novel  approaches  towards  elu- 
cidating the  physics  of  these  devices  and  generating  efficient  techniques  for 
calculating  the  eigenmodes  and  eigenvalues  of  unstable  optical  resonators. 
Below,  we  have  sunmarized  the  results  and  related  our  work  to  that  of  others 
in  the  field.  We  also  discuss  the  relevance  of  our  approach  to  problems 
involving  unstable  optical  resonators  with  gain  media. 


To  clarify  the  motivation  for  our  research,  we  first  briefly  review 
1 • 2 

the  work  of  P.  Horwitz  and  R.  Butts  and  P.  Avizonis  . These  workers  have 
developed  asymptotic  approaches  to  unstable  resonators  which  we  have  been  able 
to  complement  and  extend. 

In  1973,  Horwitz  presented  a very  interesting  asymptotic  approach  for 
treating  the  electrodynamics  of  large  Fresnel  number  resonators  with  rectangu- 
lar geometries.  The  principle  result  of  this  calculation  was  the  discovery  of 
a set  of  rapidly  varying  diffraction  functions  which  could  be  used  to  obtain 
the  empty  resonator  eigenmodes.  Specifically,  the  resonator  eigenmode  was 
written  as  the  sum  of  a fundamental  spherical  wave  plus  a series  of  edge  dif- 
fracted waves.  The  resultant  series  converged  rapidly  for  large  magnifica- 
tions (M)  with  the  convergence  becoming  poorer  as  M decreased  towards  unity. 
Recently,  Butts  and  Avizonis  extended  Horwitz's  analysis  to  resonators  with 
circular  mirrors.  The  eigenmode  was  once  again  written  as  a spherical  wave 
plus  a series  of  edge  diffracted  waves  appropriate  to  cylindrical  coordinates. 

Despite  the  very  interesting  results  of  these  authors  a number  of 
difficulties  still  remained.  In  particular, 

(1)  Convergence  was  best  for  large  M.  In  their  analysis  Butts  and 
Avizonis  found  that  this  approach  failed  entirely  for  M < 1.3 
with  N * 10. 
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(2)  For  both  cylindrical  and  rectangular  geometries  the  edge  dif- 
fracted waves  contained  unphysical  infinities  at  the  shadow 
boundary  of  the  resonator. 

(3)  It  is  difficult  to  extend  this  approach  to  problems  involving 
driven  laser  resonators,  especially  in  the  non-linear  regime. 

It  is  clear  that  an  asymptotic  approach  to  solving  the  resonator  equ- 
ation is  relevant  in  the  limit  of  Fresnel  numbers  where  numerical  solutions 
are  difficult  to  obtain.  One  would  like  to  eliminate  the  unphysical  singulari- 
ties that  appear  in  the  existing  asymptotic  solutions,  and  extend  the  applica- 
bility of  asymptotic  solutions  to  devices  with  small  magnifications.  Finally, 
it  is  necessary  to  understand  the  physical  consequences  of  the  approximations 
involved  in  the  asymptotic  approach  to  adequately  treat  the  more  complicated 
case  of  a resonator  with  an  active  gain  medium. 

Under  AFOSR  support,  we  have  developed  an  asymptotic  [and  analyticj 
approach  to  unstable  resonators  that  does  not  suffer  from  the  above  diffi- 
culties. In  particular, 

(1)  Our  approach  demands  only  that  the  Fresnel  number  of  the  reson- 
ator be  large,  is  valid  for  all  magnifications,  and  is  not  de- 
pendent on  the  specific  symmetry  of  the  systems.  Approximations 
are  made  on  the  integral  equation  itself,  rather  than  the  basis 
functions,  leading  to  results  which  are  generally  applicable. 
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For  a fixed  N,  the  solutions  to  the  resonator  equation  divide 
into  two  classes.  For  curved  mirrors  with  M>l+l//fT  the  modes 
are  diffraction  dominated,  and  our  work  reduces  to  that  of  Butts 
and  Avizonis  (for  circular  mirrors)  and  Horwitz  (for  rectangular 
devices) . 

For  nearly  flat  mirrors  with  M<1+1//?T,  the  solutions  are  domin- 
ated by  the  waveguide  nature  of  the  cavities,  and  are  given  by  a 

3 

smooth  core  term  analogous  to  that  obtained  by  Weinstein  , 
together  with  a small,  rapidly  oscillating  function  representing 
the  effects  of  edge  diffraction.  The  core  term  contains  the 
diffraction  loss  effects,  and  the  oscillating  term  averages 
essentially  to  zero. 

In  our  approach  to  resonator  electrodynamics  the  fluctuation 

effects  of  diffraction  appear  in  the  form  of  edge  waves.  As 
discussed  in  Appendix  I as  well  as  references  (4)  and  (5),  each 
edge  wave  is  the  Fresnel  diffraction  pattern  of  a plane  wave 
around  the  mirror  edge  which  has  propagated  k times  across  the 
resonator.  Furthermore,  the  edge  waves  used  by  Horwitz  (for 
strip  systems)  and  Butts  and  Avizonis  (for  cylindrical  reson- 
ators) are  the  leading  terms  in  the  asymptotic  expansions  of  the 
Fresnel  diffraction  functions. 

No  unphysical  infinities  appeared  in  our  treatment  of  the  reson- 
ator eigenmodes,  because  the  diffraction  functions  are  given 


exactly  in  terms  of  Lommel  functions.  The  leading  terms  of  the 
asymptotic  expansions  of  the  Fresnel  diffraction  patterns  con- 
tain singularities  on  the  shadow  boundary. 

(6)  As  a critical  test  of  our  approach  we  have  applied  it  to  a mar- 
ginally stable  (M=l ) cylindrical  resonator  with  a Fresnel  number 
of  ten.  The  resultant  eigenvalues  and  eigenmodes  were  in  excel- 
lent agreement  with  computer  generated  solutions  of  the  Fox-Li 

6 7 

resonator  equation.  ’ This  work  was  reported  in  references  (4) 
and  (5)  which  are  enclosed  with  this  report. 

(7)  For  the  waveguide  region  M<l+1//N,  we  have  demonstrated  that  in 
the  limit  the  core  term  reduces  to  the  closed  resonator  or 
waveguide  solution,  and  the  edge  diffraction  term  vanishes.  We 
emphasize  that  previous  asymptotic  theories  do  not  apply  to  this 
regime,  and  for  large  N reduce  to  geometric  optics  rather  than 
closed  resonator  solutions. 

(8)  For  the  diffraction  dominated  region,  M>l+I/vTTi  the  core  term 
reduces  in  the  limit  of  large  N to  the  fundamental  spherical 
wave  used  by  Butts  and  Avizonis  (for  circular  systems)  and 
Horwitz  (for  rectangular  systems). 

% 

In  Appendix  I we  have  outlined  the  most  salient  features. of  our  treat- 
ment of  empty  resonators  in  a form  that  is  applicable  to  all  values  of  M.  The 


« 


8 


basic  elements  of  the  theory  are  surrmari zed  in  Table  (1-1)  for  cylindrical  and 
rectangular  geometry. 

The  essence  of  our  approach  is  to  approximate  the  Green's  function 
for  the  differential  equation  governing  resonator  dynamics,  and  then  to  solve 
the  approximate  integral  equation  exactly.  Consequently,  the  application  to 
gain  media  is  straightforward  since  the  gain  is  a source  term  in  the  Green's 
function  solution.  In  Appendix  II  we  outline  the  derivation  of  the  unstable 
resonator  equation  for  driven  resonators.  The  same  type  of  approximations 
made  for  empty  resonators  may  be  applied.  Approximate  solutions  both  for 
linear  gain  and  non-linear,  saturable  gain  resonators  are  also  discussed  in 
Appendix  II. 

Specific  points  of  interest  which  have  arisen  from  our  treatment  of 
driven  resonators  are: 

(1)  In  the  paraxial  approximation.  Maxwell's  equations  reduce  to  a 
form  similar  to  a diffusion  equation  in  the  transverse  direction,  with  the 
longitudinal  coordinate  z in  place  of  time.  The  Green's  function  spreads 
transversely  and  decreases  in  amplitude  as  z is  increased,  in  complete  analogy 
with  diffusion. 

(2)  The  integral  equation  for  driven  resonators  implicitly  contains 
diffraction  off  the  edge  of  a non-uniform  gain  media  in  the  same  manner  that 
the  empty  resonator  equation  describes  diffraction  around  the  mirror  edge. 
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Also,  we  expect  the  equation  to  describe  the  phenomena  of  refraction  and 
internal  reflection  at  the  boundary  of  the  gain  medium  and  empty  space. 

(3)  In  the  waveguide  region  (M<l+l/»/N)  with  a non-linear  gain 

medium,  the  resonator  saturates  very  quickly  as  single  pass  unsaturated  gain 
exceeds  diffraction  losses.  The  reason  for  this  is  that  each  ray  makes  a 

large  number  of  passes  through  the  gain  medium  before  escaping,  leading  to 
large  amplifications  and  saturation. 

(4)  For  the  diffraction  dominated  region  with  M > 1+1/.4T,  each  ray 
makes  a relatively  few  number  of  transits  before  leaving  the  cavity,  and  satur- 
ation is  expected  to  be  of  less  importance.  Different  methods  must  be  used 
for  this  region. 

(5)  We  have  shown  in  Appendix  II  how  a homogeneous  linear  gain 

medium  leads  to  the  same  integral  equation  as  that  for  an  empty  resonator, 
with  the  eigenvalue  jpp  modified  by  the  single  pass  unsaturated  gain,  x = 
exLx0iFp.  The  modes  have  exactly  the  same  form  as  those  of  an  empty  resonator. 

(6)  For  a saturable  gain  medium  in  a cavity  with  flat  mirrors  the 
resonator  saturates  for  single  pass  gain  exceeding  diffraction  losses  by  as 
little  at  .1%.  The  modes  have  nearly  the  same  form  as  for  the  empty  cavity. 

In  sumnary,  we  feel  that  our  research  illuminates  the  fundamental 

physics  that  govern  the  electrodynamics  of  unstable  resonators.  Our  methods 

are  applicable  to  different  mirror  geometries,  are  valid  in  the  waveguide 
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dominated  region  where  even  ray  optical  solutions  become  cumbersome,  and 
describe  the  resonator  with  a gain-filled  cavity.  Pi.sently  we  are  obtaining 
numerical  solutions  for  empty  resonators  with  l^l+l/ZN-.  These  solutions 
and  the  methods  used  to  obtain  them  will  be  used  to  predict  the  amplitude  and 
distribution  of  radiation  in  saturable  gain  unstable  resonators. 
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APPENDIX  I 


OUTLINE  OF  ASYMPTOTIC  APPROACHES  TO  EMPTY  UNSTABLE  RESONATORS 


As  a specific  example,  we  consider  a cyl indrical  ly  symmetric  unstable 
resonator  with  magnification  M and  a Fresnel  number  N.  The  Fox-Li  resonator 
equation  can  be  cast  in  the  following  form 


V *nl(r> 


O • l+l M 
2m  N 


{ 


(I’D 


where  the  eigenfunction  of  the  (n,Jt)  Fox-Li  transverse  mode  is  given  by  f = 

i n*» 

2 

4>nl  exp(-iirNeqr  ),  the  eigenvalue  is  given  by  Tru  = and  is  the 

Bessel  function  of  order  i.  As  a first  step  we  rewrite  the  reduced  integral 
Eq.  ( I -1 ) as 


-i*N[p2+(r/M2)]M 


\ I dPP  JJl(2irNppj)e 


(1-2) 


0 1 


We  shall  refer  to  the  first  term  on  the  right  hand  side  of  Eq.  (1-2)  as  the 
core  term  and  the  second  as  the  edge  term.  The  physical  meaning  of  these 
terms  will  become  apparent  as  we  proceed.  Next  we  make  the  following  approxi- 
mation on  the  edge  term, 


( 
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2m 


;i+l 


N \ dpp 

r 


(2*Np£)e 


-i"N[P2+(-£)] 


»n*(p)  = V U)  F*(p.N),  ( I —3a) 


where 


F£(r,N) 


O • i+l M 

2m  N 


(2*Np^)e 


-iirN[p2+(^) 


(I'3b) 


There  are  three  features  that  should  be  noted: 

(1)  Using  asymptotic  techniques,  one  can  demonstrate  that  equation 

(I-3a)  is  valid  to  order  (N)-1  « 1.  For  r<l,  the  integral  Eq. 
(I-3a)  is  dominated  by  the  contribution  from  the  lower  endpoint 
of  the  contour,  and  <i>np(p)  may  be  replaced  by  ♦np(l)  to  lowest 
order  in  N-'*'. 

(2)  Ft(r,N)  is  the  *th  azimuthal  Fresnel  diffraction  pattern  asso- 

ciated with  the  mirror  edge.  For  4=0,  it  is  the  compliment  of 
the  s-wave  diffraction  pattern  of  a circular  hole,  which  was 

discussed  extensively  by  Lofnnel . 

(3)  The  integral  in  Eq.  ( I -3b)  can  be  analytically  expressed  in 

terms  of  a rapidly  convergent  Neumann  series  of  Bessel  functions. 
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The  core  term  in  Eq.  (1-2)  is  dominated  by  the  stationary  phase  point  for  r<M, 
and  is  of  order  <}>np  (p/M).  The  edge  term  which  is  of  order  \//  N arises  from 
diffraction  and  varies  rapidly  in  space.  To  obtain  a physical  understanding 
of  the  core  term,  we  neglect  diffraction  entirely  and  solve  equation  (1-2) 
exactly.  Following  this,  we  shall  incorporate  diffraction  and  solve  equation 
(1-2)  using  the  approximation  (I -3a).  The  resultant  eigenmodes  and  eigen- 
values are  exact  to  order  (1/N). 

For  cylindrically  symmetric  systems,  we  find  in  the  absence  of  dif- 
fraction 


where  r is  the  gamma  function,  ^F^  is  the  confluent  hypergeometric  function, 
and  s can  be  complex. 

There  are  several  features  of  Eq.  (1-4)  that  are  of  interest: 

(1)  The  eigenvalue  has  the  geometric  optics  value  of 
unity  for  s=0  and  the  loss  per  bounce  is  1 — L 
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For  M=l,  b=0  and  Eq.  (1-7)  reduces  to  a single  Bessel  function 
of  complex  argument, 

4*at(r)  = Jt(ar)  . (1-9) 

(4)  For  a strip  resonator  one  obtains  parabolic  cylindrical  func- 
tions. For  M^l,  the  asymptotic  expansion  of  these  modes  are 
the  geometric  optics  solutions  which  correspond  to  the  funda- 
mental spherical  wave  used  by  Horwitz  in  reference  (1).  For 
M=l,  they  reduce  to  the  modes  of  a closed  strip  resonator, 
sin(ax) . 

If  diffraction  is  included,  then  to  order  (1/N ) the  empty  resonator 
electrodynamics  are  specified  by 


»+1  r r "ilrN  ^p2+^  3 

2wil+1N^j  dpp  Jt(2irNojj)e  M ♦oji(p) 


The  exact  solution  of  Eq.  (I -10)  is 


ni 


w IS 

(r)  = ♦ pjU)  (it-)  fj,  * fT-  ) 

nl  k=l  ni  1 Mk  Vl 


♦nt(l)F(r,N). 

(1-10) 


d-11) 
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where 

k-1 

Mk-1  E £ 72J  (N12) 

j=0 

and  the  \nl  are  obtained  from  continuity  of  the  eigenmode  at  the  mirror  edge. 
Eq.  ( I -1 1 ) corresponds  to  a coherent  superposition  of  diffracted  waves  which 
propagate  back  and  forth  across  the  resonator  and  are  modified  by  the  magni- 
fication. The  usefulness  of  this  solution  is  set  by  the  convergence  of  the 
series  which  is  turn  is  determined  by  the  resonator's  magnification.  There 
are  in  fact  several  different  regions  of  interest.  These  are:  (1)  M » 1, 
(2)  M=1  and  (3)  l<M<l+l//fT. 

Case  (I)  M»1 

For  this  case  convergence  is  excellent.  To  use  equation  (I -11)  we 

r*  m 

note  that  for  sufficiently  large  k,  F (— r-,  ) is  independent  of  both  r and 

M Mk-1 

as  k,  i.e., 

^i^T’  U — )^=  no.^pr  ) if  k > k0  (1-13) 

1 Mk  nk-l  \ -1  0 


f 
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and  one  can  then  rewrite  Eq.  (I -11)  as 


F • N_  ) + 

' fl  V 1/  9 M l 

1 Mk  n k-1 


♦ ..(1) 


ni 


F,(0;- 


V1 


) S (f- ) 


V1 


ni 


k k 

• •„,»)  £ (f->  ft-  > + <xm'l) 

k=l  m M k-1  kQ-l 

(1-14) 

The  first  term  in  Eq.  (1-14)  corresponds  to  a sequence  of  edge  diffracted 
waves  and  the  second  to  a spherical  wave.  We  obtain  the  Butts  and  Avizonis 
description  if  the  edge  waves  are  given  by  their  asymptotic  expansion. 

Case  (2)  M=1 

For  this  case,  one  requires  all  of  the  terms  in  the  series  ( I -1 1 ) and 
the  techniques  of  reference  (2)  can  not  be  used.  Eq.  (I— 11 ) becomes 


ni 


k 

(r)  = 0nl(l)  £ (t5-) 


x i F (r;N/k) 
k=l  nt 


(1-15) 
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which  consists  of  a coherent  superposition  of  Fresnel  diffraction  patterns 
that  propagate  1,2, ...,k,...  times  across  the  resonator.  To  obtain  solutions 
from  Eq.  (1-15)  we  note  that  diffraction  will  perturb  the  low  lying  transverse 
modes  only  weakly  from  the  closed  resonator  form  [e.g.,  the  loss  per  bounce  is 
only  a few  percent  for  the  TEMg  Q mode].  Thus,  we  assume  that  all  but  the 
first  term  in  Eq.  (1-15)  coherently  sum  to  form  a core  term  which  is  conveni- 
ently represented  in  terms  of  a Fourier-Bessel  series,  t 

V(r>  - % cjr’o.Cj,  r)  * -±-  F((riN),nt(l),  (1-16) 

j-i  ni 

where  the  a.  are  determined  by  requiring  that  the  Bessel  functions  be  ortho- 
J * 

gonal  within  the  resonator  and  the  C^  are  obtained  by  setting  Eq.  (1-16) 
equal  to  Eq.  (1-15).  Convergence  is  rapid  and  only  a few  Bessel  functions  are 
required.  This  particular  case  is  discussed  in  references  (4)  and  (5).  Note 
that  the  eigenmodes  consist  of  a slowly  varying  core  term  with  small  oscilla- 
tions from  the  Fresnel  diffraction  pattern  of  the  mirror  edge.  The  core  term 
can  be  shown  in  the  limit  of  large  N to  approach  a single  Bessel  function 

J.(b„„r),  where  b . is  the  nth  zero  of  the  ith  Bessel  function. 
i nt  ' nt 

Case  (3)  1<M<1  + — 

/N 

A spatially  varying  core  term  which  is  not  described  by  geometric 
optics  remains  and  a number  of  diffraction  functions  are  required.  Recalling 
that  the  solutions  Eq.  (1.4)  for  M^l  consist  of  confluent  hypergeometric 
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functions  one  is  tempted  to  use  a series  of  such  functions  for  the  core  term. 
For  computational  purposes,  it  is  easier  to  work  with  Bessel  functions  and  one 
can  express  the  hypergeometric  functions  in  terms  of  Bessel  functions  via  a 
Neumann  series  as  in  (1-7).  One  then  sets 


'ni, 


(r) 


P 

= 2 
m=l 


2 

j=o 


'j* 


J*+j(ar) 


ko 

£ 

k=l 


(t1-) 


ni 


F‘(? 


nr-) 


■k-i 


and  proceeds  from  there.  Numerical  calculations  indicate  that  as  M increases, 
the  core  term  declines  in  magnitude  and  moves  outwards  towards  the  shadow 
boundary.  This  case  is  presently  under  study. 

One  can  carry  out  an  entirely  analogous  procedure  for  strip  reson- 
ators. For  this  case,  the  Fresnel  diffraction  patterns  can  be  expressed  in 
terms  of  the  well  known  sin  and  cosine  integrals.  We  find  the  analysis  to  be 
entirely  analogous  to  the  cylindrical  case.  The  results  are  summarized  below 


APPENDIX  II 


UNSTABLE  RESONATOR  EQUATION  WITH  GAIN 


The  Unstable  Resonator  Equation  with  or  without  gain  is  derived  from 
the  Maxwell  equations  with  the  boundary  condition  that  the  transverse  field 
vanish  on  the  mirror  surface.  The  time  dependence  is  separated  as  follows: 

t (r,t)  -E(r)  e-1wt  . (Il-l) 


Then  the  Maxwell  equations  become,  assuming  u=u0. 


, 2%  /*►»  k ■>  ■* 
k )E(r)  = -7-  P(r) 


> 


( 1 1 -2 ) 


where  k=w/c  and  P is  the  electric  polarization.  For  a resonator,  the  trans- 
verse field  will  consist  of  a left  moving  and  a right  moving  wave: 


= £ [i£(r,z)eiie 

i=0  L 


e~ikz  + C r.L-2)eUe  eikz] 


( 1 1 -3 ) 


where  cylindrical  coordinates  appropriate  to  circular  mirrors  have  been  used, 
£ £ 

the  amplitudes  and  are  slowly  varying,  and  C is  a constant  to  be  deter- 
mined. Assuming  that  P-j-(r)  satisfies  a similar  relation,  and  making  the  par- 

O 

axial  approximation  , we  find  that  each  amplitude  in  Eq.  ( 1 1-3)  satisfies  an 
equation  of  the  form 
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(H-4) 


«i  + 77?-7-»kir)  >M(r’Z’)  = 


k2 


PL,R^r>z'^ 


where  z'  = z for  £|_  and  z'  = L-z  for  zR. 

This  equation  may  be  solved  using  Green's  functions.  The  polariza- 
tion on  the  right  is  regarded  as  a source,  and  the  solution  is  determined  by 
boundary  conditions  on  the  mirror.  The  field  on  the  mirror  is  naturally 
assumed  to  vanish  away  from  the  edge: 


:(r,z' ) =fa(r,o)  r<a 
t o r>a 


(1 1-5) 


This  condition,  together  with  the  paraxial  approximation,  are  the  only  approxi- 
mations which  are  made  in  the  derivation  of  the  resonator  equation.  It  is 


assumed  that  both  the  field  and  its  derivations  vanish  at 


Then  using 


Hankel  transforms,  or  results  for  the  diffusion  equation,  the  Green's  function 
for  Eq.  (1 1-4)  is  found  to  be 


-ik(r2+r'2) 

G(r,z;r'z')  -fl£e  2‘z-z'> 


J (krr'/z-z‘)u(z-z') 


(1 1-6) 
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where 


u(x)  = 1 for  x > 0 (I  I -7 ) 

0 for  x < 0 


Treating  the  polarization  as  a source,  and  using  the  boundary  condi- 
tions given  by  Eq.  ( 1 1 -5 ) , the  solution  to  Eq.  (1 1-4)  is  given  by 


'L.R 


_ k ,4+1  Q-ikr  /2z' 


(r,z‘ ) = ^ i‘+i  e 


■f 


pdpe' 


-ikP^/2z* 


Jt(krp/z'  )t[>r(p»0) 


00 

{ pdp e-mrVl/ZU'-z' ' )j  (krp/(z..z.  ■ ))e\,  >z. . , 

J L,R 


(11-8) 

The  plane  z'=0  coincides  with  the  mirror  surface  only  for  flat 
mirrors.  Assuming  a curvature  R for  the  mirror,  and  R»a,  t*  R(r,o)  is  re- 
lated to  an  amplitude  ♦ on  the  mirror  surface  by 


(r,o)  = e 


ikr2/2R 


L,R 


(r) 


(11-9) 


This  assumption  may  be  seen  to  yield  the  correct  results  for  curved  mirrors  by 
obtaining  the  Green's  function  in  spheroidal  coordinates.  The  Green's  func- 
tion in  spheroidal  coordinates  has  the  same  form  as  Eq.  ( 1 1-6 ) . In  the  limit 


of  slightly  curved  mirror,  R/L»l,  the  following  results  are  obtained.  Substi- 
tuting Eq.  (II -9 ) into  Eq.  (II -8 ) , we  arrive  at  the  amplitude  at  the  point 
r,z‘  given  the  amplitude  on  the  mirror  surface  and  the  polarization: 


lR(r’z,) 


a 

k .i+l  r 
= J 


Pdpe-1k(r2+‘,2>/2zVk‘'2/2R 


Ja(krp/z*)4>LjR(p) 


z » 


^ pdp  e'ik(r2+')2)/2(z,-z,,)j)i(krp/(z1-z1‘))Pf>R(p,z**) 


(11-10) 

To  detennine  <t>,  we  develop  an  eigenvalue  equation.  Letting  z’=L,  we 

£ 

first  use  Eq.  ( 1 1 -1 0 ) to  determine  z (r,L),  which  is  then  related  to  the  ampli- 
tude <t>  on  the  mirror  surface  with  Eq.  ( 1 1 -9 ) . Using  the  boundary  condition 
that  the  total  field  given  by  Eq.  (1 1-3)  vanishes  on  both  mirrors,  we  obtain 


lL(r,L)e"ikL  + CiR(r,o)e+ikL  = 0 


iL(r,o)  + CtR(r,L)  = 0 


(11-11) 


Define  the  following  integral  transforms: 


a 

KL,M  f ■ T \ »d«e-1k<M2+1>t''2*»2)/4ML  Jt(krp/L)f(») 


(11-12) 
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where  the  magnification  is  related  to  mirror  curvature  by 


= 1-L/R 


(H-13) 


Using  Eqs.  ( 1 1-9)- (1 1-13) , we  arrive  at  two  coupled  equations  for  the 
resonator  with  equal  mirrors: 


XR  ^R^  ~ KL,M<Mr)  " 


ik  -ikr  /2R‘ 
2e  'o 


dz  KL-z,lPL(r’z) 


XL  Mr)  = KL,M  ' 2e„ 


7Tq  e“1kr2/2R  ) dz  \-2,lPR(r>z) 


(11-14) 


where  x^  = -Ce^ 

x,  = -1/C 


(1 1 -1 4a ) 


For  equal  mirrors,  symmetry  implies  C=(-l)(l41  e“^kL,  XR=XL=X*  4>R=<j»L=<t> , and 
PL=pR=p.  In  this  case,  the  eigenvalue  equation  for  the  amplitudes  on  the 


mirror  edge  is: 


x*  = K.c 


ik  i kr2/2R 
2eo 


dz  KL-z,1P 


(11-15) 
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where 


x-  e-f(kL-"<<>  . 


(11-16) 


Eq.  (11-15)  is  first  solved  to  yield  the  amplitude  4>  on  the  mirror,  and  the 
eigenvalue  x which  determines  the  single  transit  gain  or  loss.  The  amplitude 
at  point  z1  may  then  be  found  by  using  Eq.  (II -10 ) . Finally,  the  total  field 
at  any  point  is  found  by  adding  the  left  and  right  moving  waves  as  in  Eq. 

( I I -3 ) . Since  P is  in  general  a function  of  the  total  field,  Eq.  (11-15)  is 
usually  difficult  to  solve. 


est: 


Before  discussing  specific  examples,  we  note  three  points  of  inter- 


(1) For  P=0,  Eq.  (11-15)  reduces  to  the  usual  Fresnel-Ki rchoff  equa- 
tions for  empty  resonators. 

(2)  Diffraction  off  the  edge  of  the  gain  media  is  included  if  we 
assume  that  P=0  for  r>a.  Then  Ka  rather  than  K°“  is  used,  analogous  to  the 
case  of  mirror  diffraction  for  empty  resonators. 

(3)  The  transform  of  <t>  in  Eq.  (11-15)  depends  on  M as  appropriate 
for  a mirror  term,  while  the  transform  of  P is  the  usual  Green’s  function 
which  is  independent  of  mirror  curvature. 
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V 

We  first  consider  the  case  of  linear  gain,  described  by  P=xE  where  x 
is  a constant.  We  use  the  following  properties  of  the  integral  transforms  for 
unit  magnification: 

r 

L 

I KL-z  Kz  = KL  ' (H-17) 

We  find  that  the  eigenvalue  equation  reduces  on  interation  to: 

H = exp(-ikxL/2e0)K^<f  (11-18) 

which  is  identical  with  the  empty  resonator  equation  except  for  the  additional 
phase  factor.  Consequently  the  eigenvalues  will  differ  by  the  phase  factor, 
and  the  modes  will  be  identical.  This  is  the  usual  case  where  single  pass 
gain  is  multiplied  by  the  eigenvalue  for  empty  resonators  to  yield  the  total 
gain,  for  which  no  stable  solution  exists. 

To  qualitatively  discuss  solutions  for  saturable  gain  media,  we 
choose  a polarization  of  the  form 

-m2/Es 

P = xEe  , (H-19) 

which  leads  to  analytic  results  if  various  approximations  are  made. 
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(1)  Kj^_z  in  the  gain  integral  is  evaluated  in  the  geometric  optics 

1 imit. 

(2)  The  gain  is  assumed  to  be  large  enough  so  that  saturation 
occurs,  and  the  amplitude  for  flat  mirrors  is  nearly  that  of  an  empty  reson- 
ator. 


<t>  * E0Jj,(ax)  » (11-20 

where  a is  determined  in  refs  4 and  5 and  the  equilibrium  amplitude  EQ  is  to 
be  determined. 

(3)  The  phase  changes  induced  by  polarization  in  the  modes  and  eigen- 
values are  neglected. 

For  equilibrium  to  occur,  the  radiation  in  the  cavity  neither  grows 
or  decays  and  |X|=l.  The  eigenvalue  equation  is  then  easily  solved  for  EQ. 
Leaving  out  the  details,  the  solutions  are  presented  below  in  figure  (I I -1 ) 
where  ( E0/ Es)  is  graphed  as  a function  of  gain. 

From  the  figure,  we  see  that  the  cavity  saturates  (E0^ES)  for  single 
pass  unsaturated  gain  exceeding  diffraction  losses  by  as  little  as  .1%.  The 
eigenmodes  were  found  to  be  essentially  identical  with  the  empty  resonator 
modes,  with  the  greatest  differences  near  the  mirror  edge. 
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( 


SATURATION  PARAMETER. 
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